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Abstract: Using the technology of the CAESAR approach to resummation, we examine 
C^^ \ the jet-veto efficiency in Higgs-boson and Drell-Yan production at hadron colliders and 

show that at next-to- leading logarithmic (NLL) accuracy the resummation reduces to just 



in 

«^ ' a Sudakov form factor. Matching with NNLO calculations results in stable predictions 

^D ■ for the case of Drell-Yan production, but reveals substantial uncertainties in gluon-fusion 

Higgs production, connected in part with the poor behaviour of the perturbative series for 
the total cross section. We compare our results to those from POWHEG with and without 
reweighting by HqT, as used experimentally, and observe acceptable agreement. In an 
appendix we derive the part of the NNLL resummation corrections associated with the 
radius dependence of the jet algorithm. 
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1 Introduction 

One of the most active topics currently in particle physics is the search for the Higgs 
boson [1-3]. This scalar boson is the one particle of the standard model that remains to be 
found, and its discovery would provide the most direct evidence to date that the existence 
of a Higgs field [4-8] is responsible for electroweak symmetry breaking. 

Several production and decay channels are used to search for the Higgs boson at the 
LHC and Tevatron, with varying sensitivities across the range of Higgs-boson masses allowed 
in the standard model. One channel with sensitivity over a broad range of masses is gluon- 
fusion production followed by decay to a W~^W~ pair, one of which may be offshell. For 
a Higgs mass of about 125 GeV, as hinted at by current data [1, 2], this channel is one 
of several that is expected to be clearly observable with forthcoming data. An accurate 
understanding of its cross section will therefore be crucial in constraining the Higgs-boson 
couplings in the event of a discovery. 



Extensive work has been performed over the past decade to precisely predict the cross 
sections and branching ratios for all the main production and decay channels, as reviewed for 
example in [9, 10]. In the case of Higgs boson production via gluon fusion with subsequent 
decay to two W bosons, a feature that is of particular importance in discussing cross-section 
determinations is that it is customary for the LHC experiments to separately treat events 
according to their different jet multiplicities [11, 12]. This is because events with no jets 
are relatively free of backgrounds other than W^W~ production, while events with one jet 
have additional backgrounds from it and Drell-Yan production, and events with two jets or 
more also receive enhanced signal contributions from vector-boson fusion production. 

Of particular usefulness for the problem of evaluating cross sections in gluon-fusion 
production with different jet multiplicities, are fixed-order calculations that allow one to 
place arbitrary cuts on the final state, up to next-to-next-to-leading order (NNLO) for in- 
clusive production [13, 14] and NLO for the production in association with one [15-17] or 
two jets [18, 19]. However, the transverse momentum (pt) thresholds used for identifying 
jets are usually well below the Higgs mass (Mh)- This results in the appearance of logarith- 
mically enhanced terms at all orders of perturbation theory a" In " M^/pt, which can spoil 
the convergence of fixed-order truncations of the series. Such problems arise also when ex- 
amining other related final-state observables, such as the Higgs-boson or Z-boson transverse 
momentum, and the "beam thrust" [20], for which resummations have been performed to 
next-to-next-to-leading-logarithmic accuracy (NNLL) and combined with NNLO [21-24], 
or the transverse energy fiow, resummed to next-to-leading-logarithmic [25] accuracy. 

Perhaps surprisingly, however, to date no resummation has been performed for Higgs- 
related final-state definitions that involve jet finding. This may be because jet related 
observables are less inclusive than those resummed at high accuracy so far, and therefore 
require a more detailed understanding of how multiple emissions affect the value of the 
observable. For example in notable cases, e.g. the y23 jet resolution parameter of the 
exclusive fct-algorithm in e^e~ or hadronic collisions [26, 27], the existing NLL answer has 
yet to be expressed in anything other than a numerical form [28, 29]. Therefore, instead, 
the matched NNLL+NNLO results for other observables like the Higgs pt have been used to 
reweight NLO hadronic event generators, MC@NLO [30] and POWHEG [31] (those generators 
have been also used in standalone form). Insofar as a jet veto differs from a Kiggs-pt 
veto starting only from order a^ relative to the Born process, the reweighting procedure 
should provide reasonable predictions. On the other hand, the modelling of inclusive Higgs 
production processes in MC@NLO or POWHEG is such that the a^ difference between the jet 
veto and Higgs pt veto is not correctly included. Consequently some of the accuracy of the 
NNLL+NNLO calculation is lost in the reweighting procedure. 

The purpose of this paper is to examine the resummation, matching to fixed order and 
resulting phenomenology directly for the jet veto observable itself. Because a jet definition 
is a less inclusive observable than the Higgs-boson pt or the beam thrust, the resummation 
cannot be immediately reduced to standard forms such as [32]. However, at NLL accuracy, 
to which we shall mostly limit ourselves here, one may make use of the computer automated 
expert semi-analytical resummer (caesar) [33] to perform the resummation. The result 
turns out to be rather simple and straightforward also to understand analytically. 



For Higgs production, throughout this article, we will use the large rritop approximation. 
This does not affect the resummation at our accuracy but is relevant when we combine the 
resummation with fixed order calculations. 

In parallel with the Higgs-boson case, we will also examine the jet-veto resummation, 
matching and phenomenology for Drell-Yan production, which has been argued to provide 
a control case, even if, as we shall see, the issues that arise in the Higgs-boson case are not 
directly mirrored in Drell-Yan production. 

2 The jet-veto efficiency 

We consider the production of a Higgs or a Z boson accompanied by N extra QCD partons 

P1,...PN, 

pp —?' H +pi + . . .pN , and pp ^ Z + pi + . . .p^ ■ (2.1) 

A jet-veto condition is imposed by clustering the events into jets using a suitable hadron- 
collider jet-definition (JD) and requiring that the event has no jets with transverse mo- 
mentum above a certain threshold, typically in the range of 25 — 30 GeV. To define the 
jets, the LHC experiments usually use the anti-/cj algorithm [34], which repeatedly merges 
the pair of particles with smallest distance measure dij = min(p^ ,pZ )Ai??/i?^, unless 
there exists a particle with a dis = p^ value that is smaller, in which case i becomes a jet 
and is removed from the list of particles. Here ARf- = (j/j — j/j)^ -|- {(j)i — (pj)'^, where yi 
and (j)i are respectively the rapidity and azimuth of particle i. The parameter R sets the 
typical angular reach of the jet definition and is often referred to as the jet radius. After 
the clustering, one may choose to consider all jets, or alternatively only those within some 
limited rapidity range, in reflection of actual experimental acceptances. 
The cross section with a jet- veto is defined as 

/J 
'^^^d^ 0(Pt,veto - Ptji (Pi' • • • ^Pn)) , (2-2) 

where ptji is the transverse momentum of the hardest (highest pt) of the jets found in the 
event, da^ denotes the partonic cross-section to produce a Higgs or a Z boson accompanied 
by N extra partons and d^N is the corresponding phase space. 
It is also useful to consider the jet-veto efficiency defined as 

e(Pt,veto) = , (2.3) 

where a is the total cross-section. The veto efficiency is of interest because it essentially 
encodes just the information about the Sudakov suppression associated with forbidding 
radiation of jets. In contrast, the vetoed cross section mixes in also the physics that deter- 
mines the total cross section. Thus, in the absence of a veto, pt,vcto = oo, the efficiency is 
exactly 1 and one can reliably discuss small departures from e = 1 as pt,veto is reduced. In 
the vetoed cross section, T,{pt^vcto), it is harder to disentangle those effects from uncertain- 
ties on the total cross section. Later, we will argue that even for e substantially below 1, 



it makes sense to treat the efficiency and total cross section as independent quantities, and 
that the uncertainties that govern them are relatively uncorrelated. 

Each of S(pt,veto)) cr and e(pt,vcto) has a fixed-order perturbative expansion, which we 
write as 

S(pt,vcto) = So(pt,voto) + Si(pt,vcto) + S2(pt,voto) + • • • , (2.4a) 

e(pt,voto) = eo(Pt,vcto) + ei(pt,vcto) + e2(pt,veto) + • • • , (2-4b) 

O- = (To + (Ti + (T2 + . . . , (2.4c) 

where the index i signifies that the contribution is proportional to a* relative to the Born 
term. The properties So(pt,voto) = ctq and eg = 1 follow from the fact that no jets are 
present at Born level. 

3 NLL Resummation 

The jet-veto efficiency is nothing but the cumulative distribution of the transverse momen- 
tum of the highest-pi jet, ptji- In order to have a dimensionless observable, we choose to 
divide it by the corresponding boson mass, Mb = Mz or Mh, 

V{pi,p2;ki...kN) = ^^ — — . (3.1) 

Mb 

Here pi and p2 denote the incoming partons entering the hard scattering (after any initial 
state emission), ki . . . k^ are the momenta of the final-state QCD partons. The momentum 
of the boson can be obtained using energy-momentum conservation. 

We present here a next-to-leading logarithmic (NLL) resummation of the jet-veto ef- 
ficiency, i.e. we resum all logarithms in lne(pt,vcto) up to a^L"^, with L = hi{M b / Pt.wcto) ■ 
This observable is within the scope of the computer automated resummation program CAE- 
SAR [33]. CAESAR is a program that, given a computer subroutine for an observable, 
automatically performs a numerical analysis of its behaviour with respect to multiple soft 
and collinear emissions. From this analysis it establishes whether the observable belongs 
to a broad class for which it is able to perform NLL resummations, and if so expresses the 
result for the resummation in terms of a "master formula" together with various numerically 
determined input parameters. The physics framework that underlies CAESAR can also be 
used in an analytic context as we shall do, briefiy, here. 

The first element of the analysis of an observable is to establish its dependence on the 
kinematics of a single soft and collinear emission. For the case of the jet veto, we simply 
have 

n{p},A.) = A, (3.2) 

where kt is the transverse momentum of the emission with respect to the beam. This 
expression holds exactly in all (hard, soft, collinear) limits. 

The next element is to establish whether the observable is continuously global [35, 
36]: essentially, the observable should be non-zero for any emission with finite energy and 
angle with respect to the beam; furthermore, the power-law of the observable's dependence 



on the emission energy should be independent of the emission angle. Strictly speaking, 
the first part of this condition holds only if jets are measured at all rapidities, whereas 
real experimental measurements have finite acceptance. We shall however for now work 
assuming full acceptance and return to the question in section 7. Given the first part of 
the condition, the second part holds trivially for the jet veto, and this ensures that the 
resummation is free of "non-global" logarithms and related terms. 

A second condition was dubbed recursive infrared and colUnear (rlRC) safety [33], and 
concerns the observable's behaviour in the presence of multiple emissions. Essentially it 
involves two requirements: (a) if one scales all emissions in some appropriately uniform 
manner towards the soft and collinear limit, then the observable should scale in the same 
manner; and (b) if one emission is soft relative to the others, and kept of fixed relative 
softness while scaling all the emissions, then in the scaling limit, the observable should 
not change if the soft emission is removed. Again, this condition holds trivially for the 
jet veto, as long as one uses a standard "inclusive" longitudinally invariant infrared and 
collinear-safe hadron collider jet algorithm (e.g. kt [27, 37], Cambridge/Aachen [38, 39], 
anti-fet [34]). It ensures that double-logarithmic terms exponentiate and that one can use 
an independent-emission type approximation in evaluating NLL terms. 

The master formula in the CAESAR approach then tells us that the NLL resummed 
jet-vetoed cross section is given by 

SNLL(Pt,veto) = / dxidX2 /(xi,Pt,veto) f{x2,Pt,veto) 6{xiX2S - M|) \Mb\'^ 

■ T{R'B)e-^^^P''-^'°\ (3.3a) 

where [A^^p is the full tree-level matrix element squared to produce a Z or a Higgs boson 
including the flux, and spin- and colour-averages; e~^ is a Sudakov form factor, with Rb 
a double-logarithmic function given respectively for the Z and Higgs-boson cases by 

Rz (Pt,vcto) = 2C, J^^^^^^^ ^^ ^In ^ - 4 J , (3.3b) 

Rh (Pt,veto) = 2Ca j^^^^^^^ -^^ ^In — ^^^ J , (3.3c) 

with a^^^(/c() the coupling constant in the CMW scheme [40], evaluated at scale kf. 
The quantity R'^ is equal to dRB/dln{MB/pt,vcto) and to NLL accuracy it is simply 
4CQs(pt,veto)/7rln(MB/pt,veto)- It enters in the function T{R'b), which provides the single- 
logarithmic contributions associated with the observable's dependence on multiple emis- 
sions. As derived in [33] (and references therein), T{R'b) is given by 

oo I m 2 „, roc J/- fl r2TT 



-('^.)--/"^E;^ nEfff/<.r|^ 



m=0 \i=l^,=l 



x9(l-lim^(^^^'"^(^^")'---'""-(^-"^^'), (3.4) 



where the Ki{C,iv) are parametrised momenta with azimuthal angle (j)i, rapidity yi = {—lY^£,i In 1/v 
and kti = CivMs, with the hard incoming momenta {p} adjusted to ensure momentum con- 
servation. In evaluating J-", a key point is that in the limit ?; — >• all emissions become 
widely separated in rapidity. Consequently the jet algorithm, independently of the jet ra- 
dius R and whether it's of the SISCone [41] or generalised- /cf [27, 34, 37-39] type, clusters 
each particle into a separate jet. The leading jet will then be given by the hardest of all 
the emissions, so that 

v—>0 V 

It is then straightforward to show that T{R'q) = 1. This can be verified also numerically 
with CAESAR. The result J-'(R'^) = 1 has been found before also for the resummation of 
the jet-resolution parameter in the e'^e~ Cambridge algorithm [28]. 

The feature of the jet veto observable that causes T{R'^) = 1 has consequences for 
the resummation beyond NLL. In particular certain key observable-dependent pieces of the 
NNLL resummation were outlined in the appendices of [33]. These are discussed further 
in appendix A. For the remainder of the main part of the article, however, we restrict our 
attention to NLL resummation. 

In what follows below, rather than using Eq. (3.3) directly as written, we shall reduce 
it to a form where the exponent involves just LL and NLL terms, expressed in terms of a 
logarithm L = InQ/pt.vctoj with Q the hard scale with respect to which the resummation 
is defined, as well as the MS coupling ^^(/ijj) at a given renormalisation scale ^n and PDFs 
with a hard factorisation scale ^p- It is to be understood that Q, ^_r and ^p will all be 
taken of order Mb- Varying them in the neighbourhood of Mb will allow us to probe 
the impact of subleading terms. In practice, this is accomplished by expressing the CMW 
coupling in terms of the MS coupling 

„CMW_,(,,^.^), K.C.(f4)-5„„ ,3,.) 

and using for as{kt) the following expression 

/,N oJjiR)/^ , /3ihi(l-2p)\ _ N« 1 ^fl /Qv\ 

Qs(fct) = i_2 U-"4/^i?)-g Y^ ]' P = as{nR)l3oln — , (3.7) 

with 

IICa - 2nf 17Cl - 5CAnf - SCpUf 

^° = 12^^ ^' = 24^^^ • ^^-^^ 

A closed result for the NLL resummed result can then be obtained by substituting Eqs. (3.6,3.7) 

into Eq. (3.3): 

SNLL(Pt,veto) = / dxidX2f{xi,HFe~^) f {x2, flpe'^) S{xiX2S - M%) \Mb? 

■ exp [Lgi{as{nR)L) + g2{as{f^R)L)] , (3.9) 



with 
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+ ln 



2/3o^ 



Mb 
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(3.10a) 
(3.10b) 
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Q 1-2A 



where, for Higgs production C = Ca and B = — (llC^i — 2nj)/(12Cyi), whilst for Drell-Yan 
production, C = Cp and B = —3/4, and A = ols{^ii)PqL. This form has the property 
that in the exponent there are only LL and NLL terms, whereas Eq. (3.3) effectively also 
contained some higher-order contributions. Note also that the use in Eq. (3.9) of fipe~ 
rather than pt,veto in (3-3) corresponds to a NNLL difference. 

One advantage of formulating the resummed result in terms of the variable L is that it 
is then straightforward to adapt the resummation so as to ensure a sensible behaviour even 
for pt,vcto ^ Mb- This is done by "modifying" the logarithm. 



L ^ L = - In 
P 



^Pt,vcto 



Q 



»,Max 
, Ft, veto 



+ 1 



(3.11) 



veto ^^ rt.veto 



where Ptwvto i^ the largest physically accessible jet momentum. For Q <^ pt, 
L ~ {Q/pt,vcto)^/p and thus p can be chosen so as to ensure that resummation effects vanish 
quickly above the scale Q. In practice we choose p = 5. This should help limit artifacts 
that could arise with the matching procedure introduced later on in this article. 

The above way of formulating the resummation follows in the footsteps of [42], and 
includes refinements proposed in [43]. One of these refinements is that we allow Q ^ 
Mb, which accounts for the O (1) additive freedom in defining L and the point where the 
resummation sets in. It plays an important role in uncertainty estimates. 

For completeness we explicitly write the final form we use for the resummation cross 
section as S, which is simply Eq. (3.9) with L replaced by L: 



SNLL(Pt,veto 



) = / dxidx2 /(a 



dxidx2 f [xi, jipe ^) f{x2,fJ,Fe ^) 6{xiX2S - Mb)\Mb 



exp 



Lgi{as{fiR)L) + g2ias{fiR)L) . (3.12) 



Figure 1 illustrates the resummed results for the Higgs-boson and Z-production cases. 
It shows the efficiency defined as 



eNLL(Pi 



t, veto J 



SNLL(Pt, 



veto J 



SNLL(Pt^ 



Max A 
veto/ 



(3.13) 



The vertical dashed lines mark the values of the jet-veto thresholds used by ATLAS and 
CMS experiments, of 25 GeV and 30 GeV, respectively. The thick solid (blue) line corre- 
sponds to our central scale choice jir = iip = Q = Mb/'^.. The finely hashed (blue) band 
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Figure 1. NLL resummed result for the jet- veto efficiency in Higgs (left) and Z-boson (right) 
production. The predictions are for pp collisions at a centre of mass energy of 7 TeV and use 
1VISTW2008 NNLO PDFs [44]. For Higgs production we use the large mtop approximation. 

around it corresponds to the envelope of results obtained when examining an independent 
variation of the renormalisation and factorisation scales with a constraint on their maximal 
ratio: 



{Mb Mb ,. 



/Ms Mb ,. 
MF = |— '— '^-B 



i<^ <2, 

2 ^F 



(3.14) 



The resulting uncertainty is relatively modest. More significant is the uncertainty associated 
with the variation of Q, shown by the widely-hashed (red) band, corresponding to the 
envelope of the choices Q = {Mb /4:, Mb /2, Mb} (all with fi^ = ^p = Mb/2). This 
variation probes the impact of unknown NNLL terms. In the case of Higgs production, 
one sees that the band even goes almost 10% above 1. This kind of behaviour is common 
in NLL resummation and can be traced back to the B term in Eq. (3.10b) or equivalently 
the second, negative term in the round brackets of Eq. (3.3c). The problem is less severe 
in the Z-production case, because of a cancellation due to PDF effects. These kinds of 
artifacts are one of the reasons why it is important to "match" resummations with fixed 
order predictions, especially in the region of moderate to large pt vetoes. In the next section 
we will therefore consider the structure of the fixed-order cross sections, and subsequently 
proceed to introduce our matching prescriptions and to examine the results. 



4 Jet-veto at fixed order 

The state-of-the-art of fixed-order predictions for fully differential differential partonic 
Higgs-boson and Z-boson cross sections is NNLO, i.e. the calculation of S2(pt,vcto) and 
(72, with tools like FEHIP [13] and HNNLO [14] for Higgs productions, and FEWZ ]45] and 
DYNNLO [46] for Z production. For the purpose of determining the jet-veto cross section, it 
is however also possible (and sometimes numerically cheaper) to compute only a2 with these 
NNLO tools (or from the inclusive results [47-50]), and obtain Si(pt^vcto) and S2(pt,vcto) 



from the relation 



Si(pt,veto) =(Ti + Si(pt,veto), ^i{Pt,veto) = " / ^Pt ! ■ (4.1) 

The differential distributions d'Ei/dpt and dT,2/dpt can be computed from the boson+jet 
cross sections at LO and NLO respectively, e.g. using MCFM [18, 19, 51]. We recall that, 
throughout, we use the large mtop approximation for Higgs production. 

4.1 Prescriptions for the efficiency 

There is little ambiguity in the definition of the fixed order results for the total and jet- 
vetoed cross-sections, with the only freedom being, as usual, in the choice or renormalisation 
and factorisation scale. However, given the expressions of T, and cr at a given perturbative 
order, there is some additional freedom in the way one computes the jet-veto efficiency. For 
instance, at NNLO the efficiency can be defined as 

(a)/ N ^ ^o(Pt,veto) + Si(pt,vcto) + S2(pt,veto) /. „ s 

f \Pt,vcto) — , , ) [^.za) 

but the following expressions are equally valid at NNLO, 

(b)r N _ ^o{Pt,veto) + Si(pt,veto) + S2(pt,veto) / . „, s 

e {Pt,Yeto) — ; ! l^./Dj 

O'O + 0-1 

e^ '{Pt,vcto) = IH h 2^l(^t,vcto) , (4.2c) 

since they differ relative to Eq. (4.2a) only by terms O (af), which are not under control. 
Option (a) is the most widely used, and may appear at first sight to be the most 
natural, since one keeps as many terms as possible both in the numerator and denominator. 
However, option (6) can be motivated as follows: since the zeroth order term of e(pt,vcto) is 
equal to 1, it is really only 1 — e(pt,vcto) that has a non-trivial perturbative series, given by 
the ratio of the inclusive 1-jet cross section above pt,veto, o"i_j^t (pt.vcto)) to the total cross 
section. Insofar as the 1-jet cross section is known only to NLO, in taking the ratio to the 
total cross section one can argue that one should also use NLO for the latter, i.e. 

, ^ -I '^1-jet {Pt,vcto) 
e(pt,veto) = 1 ■ • (4.3) 

It is straightforward to verify that this then leads to Eq. (4.2b). This procedure also 
coincides with the one adopted in event-shape studies in DIS and hadron-hadron collisions 
((T2 is not even known in the latter case). Option (c) is also well motivated, since it is a strict 
fixed order expansion of the ratio, so no uncontrolled terms beyond NNLO are included. 
This is the prescription that is usually adopted in e~^e~ event-shape and jet-rate studies. 

While other possibilities are also equally valid, the above three schemes capture a 
substantial part of the freedom that one has in writing the series. The size of the differences 
between them is one way to estimate the associated theoretical uncertainty and goes beyond 
the usual variation of scales. 
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Figure 2. Jet-veto efRciency for Higgs (left) and Z-boson productfon (riglit) using tliree different 
prescriptions for the NNLO expansion, see Eqs.(4.2a-4.2c). For eacli prescription, tlie thiclc solid 
line corresponds to the result obtained with /i/j = [ip = M^/z/'^, while the band shows the scale 
uncertainty as obtained with the choices of Eq. (3.14). 

4.2 Numerical results 

Figure 2 shows the NNLO results for the jet- veto efficiency in the 3 schemes discussed above. 
Each scheme is displayed as a band corresponding to the envelope of the scale variations as 
in Eq. (3.14), together with a solid line for the prediction with the central scale choice. 

In the case of Higgs production (left-hand plot) the bands barely overlap and, in the 
region of interest, pt.veto ~ 25 — 30 GeV, the three predictions differ considerably, with 
the bands spanning the range from 0.6 to 0.9 and the central values from 0.6 to 0.77 (for 
Pt,veto = 30 GeV). For yet smaller pt,veto values one sees a complete breakdown of the 
fixed-order predictions, with the efficiencies ranging from below to above 1. Curiously, 
scheme-c misbehaves also at large pt.veto ~ Mh, where it leads to an efficiency larger than 
1. 

In contrast, for Z production (right-hand plot), the situation is much more stable, with 
all three schemes overlapping and remaining sensible- looking down even to low pt,veto- 

The large differences between the Higgs and Z-boson cases call for an explanation. One 
consideration is that the gg initial state that is relevant for Higgs production implies Ca = 3 
colour factors rather than Cp = 3 colour factors. This inevitably worsens the perturbative 
convergence. The larger colour factor is however not the only culprit. Table 1 shows the 
total cross sections for Higgs and Z production calculated at different orders, as well as the 
coefficients Cj of the perturbative expansion normalised to do at the central scale choice, 
i.e. Cj = (Tj/((Toa^) such that the cross section reads 



a 



(To(l + CiOs + C2ag + . . .) 



(4.4) 



The difference between schemes a and b is in the choice of whether the 1-jet cross section is 
normalised to the NLO or NNLO total cross section. The NNLO term brings a correction 
of order 20% for Higgs production and this is reflected in the efficiencies; for Z production. 
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LO NLO NNLO ci C2 

H [pb] 5.40tl:o2 ll-96+i:95 14.7tjj 9.78 33.5 

Z [nb] 22.85^2:40 28.61?;^ 28.6lj]:^ 1-94 0.0 



Table 1. Cross-sections for Z and Higgs (large mtop approximation) production at various orders 
in perturbation theory. The central value corresponds to the default scale, fiji = fip = Mb/2^ 
Mb = Mz or Mh, the error denotes the scale variation when fin and hf are varied independently 
by a factor two around the central value, with the constraint that fiR/2 < fxp < 2/ifl, i.e. Eq. (3.14). 
Also shown are the expansion coefficients Ci of the cross-sections for the central scale fiQ = AIb/2, 
ai = aoCiaKno) . The results were obtained using the HNNLO [14] and DYNNLO [46] programs, with 
MSTW2008 NNLO PDFs and a^(Mz) = 0.11707, and a centre-of-mass energy of 7 TeV. 

LO NLO NNLO ci C2 

H [pb] 5.40+};o^ 8.99+}:?^ ^-^^^oll 5.36 10.6 

Z [nb] 22.85t|^^ 25.541°;^^ 24.5l[]:^ 0-90 "2-7 



Table 2. Same as Tab.l but with a jet-veto of 30 GeV, based on the anti-A;* algorithm with 
R=0.5. The LO cross-section is equal to that without a veto. The results have been determined by 
subtracting H + 1-jet and Z + 1-jet cross sections, obtained with MCFM, from the inclusive cross 
sections obtained with dynnlo and hnnlo, though we could equally well have used dynnlo and 
HNNLO directly. 

(T2 is essentially zero and so schemes a and b give identical results. As concerns scheme c, 
one is additionally affected by the very large size of cJi in the Higgs case [52, 53]. 

Table 2 shows corresponding cross section results with a jet-veto of 30 GeV. The Z- 
production case maintains its reasonable convergence. Interestingly, the vetoed Higgs cross 
section now also appears to have a convergence that is significantly improved with respect 
to the total cross section. As has been argued by others [54], this apparent improvement 
in convergence can be interpreted as an artifact of cancellations between two independent 
classes of large corrections: those that lead to the poor convergence of the total cross section 
and those associated with the Sudakov suppression for small jet vetoes (enhanced compared 
to the Z case due to the larger colour factors for Higgs production). The presence of such 
artificial cancellations in the cross section is one of the reasons that we prefer to consider 
the jet veto efficiency. Furthermore, as we have seen, the use of the jet veto efficiency 
provides additional handles to estimate the perturbative uncertainty through the three 
schemes discussed above. 

Given that we have indicated that the uncertainties on the Higgs jet- veto efficiency 
are associated with the poor convergence of the total Higgs cross section, some comments 
are due concerning discussions in the literature on our knowledge of the total Higgs cross 
section. On one hand there is work that aims to account for threshold and other enhanced 
terms beyond NNLO [55-58], while other work has suggested that there is a need to revisit, 
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and perhaps be more conservative in, estimates of the uncertainties on the total cross 
section [59]. Regarding improvements of the total cross section, while in general these have 
the potential to be highly valuable, a more accurate evaluation of the total cross section 
is, we believe, likely to provide little improvement for the jet-veto efficiency: indeed, our 
results indicate that there are potentially large missing O (q^) corrections to the jet veto 
efficiency and it is only through an evaluation of the NNLO H+l-jet rate that these could 
be fully constrained. 

5 Prescription for estimating total uncertainties 

From now onwards, for the purpose of quoting final efficiency estimates, we will need a 
prescription for combining different sources of uncertainty. We shall adopt the envelope 
method, used for example in [60] as well as in many other works on resummation. This 
method takes the envelope of the curves obtained from each of several different sources of 
uncertainty estimate. The logic of the scheme is that to avoid double counting uncertainties, 
at most one source should be probed at a time. Thus here at fixed order we would take the 
envelope that results from the scale variation band from one main efficiency prescription 
(scheme a) and the central values of each of the two alternative efficiency prescriptions (b 
and c). As has been found in [61], this leads to an uncertainty estimate that is quite close 
to the alternative fixed-order prescription proposed in [54], which varies the scale choice 
independently in zero, one and two-jet bins. 

In addition to estimating the uncertainty on the veto efficiency, it is also important to be 
able to estimate the uncertainty on the final vetoed cross section. The prescription that we 
propose is to treat the efficiency uncertainty as being uncorrelated with the uncertainty on 
the total cross section. This can be justified on the following grounds: the total cross section 
uncertainty is associated with our ignorance of the NNNLO corrections to its perturbative 
series. On the other hand the uncertainties in the efficiency are related to higher-order 
Sudakov type terms and to the way in which we treat the known NNLO cross-section 
contributions in calculating the efficiency. These are sufficiently different in origin that we 
believe that it is reasonable to treat them as uncorrelated. 

6 Matching NLL and NNLO results 

To obtain predictions for the jet veto efficiency that include the advantages of both re- 
summed and fixed order results, we combine them with the help of a matching procedure. 
Below, we first briefly recall the various requirements that matching should satisfy, then 
list the matching prescriptions that we adopt, and finally discuss the matched results. 

6.1 Matching prescriptions 

As is well-known there are various ways in which one can match resummed and fixed 
order calculations. Since we match the resummation to NNLO exact results, the matching 
procedure should satisfy the following three requirements: 
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1. The matched resuh should be correct up to NLL terms in the exponent and the 
expanded matched result should be correct up to and including O (a^L^""^) terms. 

2. The expanded matched result should coincide with the fixed order result up to and 
including the NNLO terms. 



3. The jet-veto efficiency should tend to one at the maximum allowed jet transverse 
momentum p^toi ^^'^ ^^^ 
should vanish at that point 



momentum p^toi ^^'^ ^^^ corresponding differential distribution de (pt, veto )/dpt, veto 



ebt,veto) = l and ^<Phl^ = q for p,^,,,, = p^-^^ . (6.1) 

"Pt,veto 

The use of the modified logarithm L, defined in Eq. (3.11) makes it relatively straight- 
forward to fulfil this condition. Note that at the LHC, ^^tveto i^ much larger than the 
Pt.vcto values that will be of interest to us. Therefore, while we take care to enforce 
Eq. (3.12), the terms involving p^to ^i^^ ^^ practice largely be irrelevant. Neverthe- 
less the use of the modified logarithm is important also because it ensures that the L 
that enters in Eq. (3.12) approaches zero for pt.vcto ^ Q- 

Even with these conditions, there is some freedom in the matching procedure. This is 
closely related to the freedom we had in expressing the efficiencies. We will therefore 
consider three matching schemes, each of which is the counterpart of one of our fixed-order 
efficiency prescriptions, so as to facilitate the comparison with the NNLO results. 
The first of our matching schemes is given by 



^matched V^t.vcto) — 
/ ^NLL(Pt,veto) \ 

V ^0 j 



O"0 + 5]i(pt,veto) + ^2{Pt,veto) — Z ( i;NLL,l(Pt,veto) + SNLL,2(Pt,veto!! 



.SNLL,l(Pt,veto) f^ , . ^+1- , 

^liPt,vetoj :: ^NLL,liPt,vcto, 



O-Q 



, (6.2) 



where SNLL,i(Pt,veto) is the term of O (a*) relative to ctq in the fixed-order expansion of 



Pt.VGtO 

„Max 
i^t.vcto 

in practice. The corresponding jet-veto efficiency is 



Snll- The factor Z = ( 1 — mIx" ) is necessary to satisfy Eq. (6.1), but is largely irrelevant 



^(a) f^ \ _ ^matched '^Pt,^ 

ax \ 
'eto/ 



AtJ-l f„ \ _ -^matched '^J^t.veto; /„ „s 

^matchedW.vetoJ - („) • [b.6) 

^matched '^^t.veto J 



It is straightforward to verify that with this matching procedure ejnatched(^t.veto) satisfies 
all three requirements listed at the beginning of this Section. 
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The second matching scheme is identical except that we replace Ti2{pt,veto) with T,2{pt 
in direct analogy with scheme b used for the fixed order efficiency definition, 



veto ; 



""matched 



(Pt.veto) 



Snll(p 



't,veto ) 



O-Q 



cro + i;i(pt,vcto) + ^2{Pt 



^5]NLL,l(Pt,vcto; / ^ , 
■ ^ I ^l(Pt,vctoJ 



with 



O-Q 



•^matched Vi^t,vcto; 



Z ( SNLL,l(Pt,veto) + SNLL,2(Pt,veto) 

. Z + l 



-Snll,i(p- 



'tjveto. 



^matched '^Pt.vcto J 
v(^) ('„Max 



, (6.4) 



(6.5) 



-'matched \Pt,\eto) 

The third scheme that we consider is formulated directly in terms of the jet- veto effi- 
ciency and is thus close in spirit to fixed-order scheme c: 



(c) 
^matched l^t,vcto 



(Pt 



(eNLL(Pt,vcto))^ 



1 + ei(Pt,veto) + e2(Pt,veto) " Z (eNLL,l(Pt,vcto) + eNLL,2(Pt,veto)) 



■Z'eNLL,l(Pt,veto) ( ei(Pt,veto) 



Z + l 



eNLL,l(P 



tjveto. 



, (6.6) 



where, in accord with earlier definitions, we have 

5^NLL(Pt,veto) 



eNLL(Pt,veto) 

eNLL,j(Pt,veto) = 

ei(pt,veto) = 

e2(Pt,vcto) = 



SNLL,i(Pt,veto) 
Si (pt, veto) 



S2(Pt,vcto) 



o-i5]i(pi 



tjveto j 



O-Q 



(7^ 



(6.7a) 
(6.7b) 
(6.7c) 
(6.7d) 



The three schemes differ in two respects: firstly they treat non-asymptotic terms (i.e. 
Pt.veto ~ Q) differently; secondly, they lead to different sets of subleading logarithms. In 
particular, schemes a and h differ at the level of N^LL terms, while scheme c differs from 
both of them by N^LL terms. We recall that NNLL differences will be probed by scale 
variations. 

The above three schemes all belong to the family of multiplicative matchings [43]. 
Other schemes exist. In particular we have also examined the \og-R scheme of Ref. [42] and 
found that its results are contained within the band defined by the three schemes above. 
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Higgs production (m^ = 125 GeV), NLL+NNLO 



Z production, NLL+NNLO 
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Figure 3. NLL+NNLO jet-veto efRciency for Higgs (left) and Z-boson productfon (right) using 
three different matching prescriptions. For each one, the thick soHd hne corresponds to the result ob- 
tained with /ijj = [ip = Q = M^jzl'^i while the band shows the scale uncertainty as obtained from 
the envelope of the choices of Eq. (3.14) and from Q-scale variation (taking Q = {AfB/4, Mb/'2., Mb} 
for /i/j = Up = M-iijzl'^,). The lower panels show the results normalised to the central scale choice 
for scheme a. 

6.2 Matched results 

Figure 3 shows the NLL+NNLO matched results for each of the three schemes discussed 
above. Each band corresponds to the envelope of renormalisation, factorisation and Q scale 
variation. Comparing to the pure resummed and fixed-order results of figures 1 and 2, one 
observes that for larger pt,vcto) the results display the features of the fixed-order results, with 
similar uncertainties and, for scheme c, the same artifact of e > 1 for certain scale choices. 
As in the NNLO case, for Higgs production the uncertainty bands from the three schemes 
mostly do not overlap, which serves to highlight the importance of the use of different 
schemes for probing uncertainties. For low pt,veto values, the results coincide with the pure 
resummation, though with uncertainty bands that are a little narrower, and that overlap 
even in the Higgs case. In the region of intermediate pt,vcto ~ 25 — 30 GeV, the central 
value from the matched calculation is closer to that of the fixed-order results than to the 
resummed results, but with a slightly reduced uncertainty, indicating that resummation is 
at the edge of its validity in this region. Here too the bands from the different matching 
schemes fail to overlap in the Higgs case. Since the bands differ at most by NNNLL terms, 
this has implications for the degree of improvement that one might expect when extending 
the resummation from NLL to NNLL accuracy. 

A direct comparison of the fixed-order and matched predictions is to be found in Fig. 4. 



15 



Higgs production (m^ = 125 GeV), NNLO v. NLL+NNLO 



Z production, NNLO v. NLL+NNLO 
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Figure 4. Comparison of fixed-order (NNLO) and matched resummed (NLL+NNLO) predictions 
for tlie jet veto efficiencies in Higgs (left) and Z production (right). The uncertainties are those 
derived from the envelope method: for both fixed order and matched results they include renor- 
malisation and factorisation scale uncertainties, as well as the scheme for defining the efhciency 
(or matching prescription). In the matched case, there is additionally the uncertainty from the 
variation of Q. The lower panels show the ratio of the results to the central matched prediction. 

Here the uncertainty envelopes encompass the full scheme a band as well as the central val- 
ues of the two other schemes. This follows the procedure outlined in section 5 and it provides 
the uncertainties that we shall use throughout the rest of the article. The efficiencies for 
the two jet- veto thresholds used by ATLAS and CMS, 25 and 30 GeV respectively, are 
summarised in table 3. For Higgs production, one observes that the absolute efficiencies 
are about 3% lower in the matched calculation as compared to the NNLO result (equiv- 
alent to a relative 5% reduction in the efficiency). The uncertainties are somewhat more 
asymmetric in the matched calculation and in particular the uncertainty towards lower effi- 
ciencies is reduced by about a factor of two. For Z-boson production, the uncertainties with 
matching are the same or larger as those of the pure NNLO result. This surprising result 
may be because the resummation explicitly involves the running of the coupling and thus, 
for low pt,vcto) directly probes the uncertainties associated with a perturbative expansion 
whose coupling constant is somewhat larger than the as{Mz/2) that appears in the NNLO 
calculation. 

7 Comparisons to other calculations 

In this section we will complement our resummed matched study so far with information 
from event generators and analytical boson-pf resummations. For brevity we concentrate 
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Higgs production {Mh = 125 GeV) 

NNLO NLL+NNLO 



Pt,veto = 25 GeV 601^1% 
Pt,veto = 30 GeV 67^^% 



57l^% 
64+^% 





Z production 


NNLO NLL+NNLO 


Pt,veto 


= 25 GeV 81+^% 81+^% 
= 30 GeV 85+}% 85+^% 



Table 3. Jet veto efficiencies and their uncertainties at NNLO and NLL+NNLO, for the values 
of pt, veto used by ATLAS and CMS, sliown for the anti-fct algorithm with R = 0.5, and based on 
MSTW2008 NNLO PDFs. 
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Higgs production (m|_| = 125 GeV), impact of hadronisation 
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Figure 5. Left: impact of a finite rapidity acceptance for jets on the jet-veto efficiency, as 
calculated with pythia 6.425. Right: impact of hadronisation and underlying event on the jet- veto 
efficiency. See text for further details. 

on the case of Higgs production, using Mu = 125 GeV throughout. 

7.1 Effects beyond the scope of matched calculations 

The matched calculation that we have performed applies to partons and assumes infinite 
detector acceptance. Experiments, however, measure hadrons, including the underlying 
event, and have limited acceptance, notably for the rapidity of the jets. 

To investigate these two effects we have taken events generated by PYTHIA 6.4 [62] 
with the Perugia 2011 tune [63]. Jet clustering for the results in this section is performed 
with Fast Jet [64]. Fig. 5 (left) shows the impact of considering jets only within some finite 
rapidity acceptance. One sees that for the choices used by ATLAS and CMS, \y\ < 4.5 and 
y < 1 5.0 1 respectively, the veto efficiencies are almost identical to those with full acceptance 
in the practically relevant range of pt.vcto- We have confirmed that this pattern holds also 
in fixed-order calculations. In contrast, if one applies a jet veto only in a more restricted 
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rapidity region, e.g. \y\ < 2.5, there are substantial differences. 

The right-hand plots of Fig. 5 show the impact of non-perturbative effects. While not 
entirely negligible, both hadronisation and the underlying event have an impact that is 
somewhat smaller than the uncertainties on the matched calculation, at least for the range 
of pt,vcto values of practical interest. This contrasts with the situation for variables that 
receive contributions from all hadrons in the event, such as the beam thrust or transverse 
energy flow, which see large contributions from the underlying event (cf. the studies in 
Ref. [25]). 

7.2 Comparisons with POWHEG and HqT 

Both ATLAS [11] and CMS [12] make use of POWHEG [65] for estimating their jet-veto 
efficiencies, interfaced to PYTHIA 6.4 (this is stated explicitly by ATLAS, we assume it to 
be the case for CMS). The ATLAS collaboration additionally reweights the events so as to 
ensure that the Higgs boson qt spectrum coincides with that of the HqT program [22] . In 
this section we compare our results to these various tools. 

We start, figure 6 (left), by comparing the HqT NNLL+NNLO result for a veto on the 
Higgs-boson transverse momentum to our NLL+NNLO result for the jet veto. Part of the 
purpose of this comparison is to examine the relative sizes of the uncertainties and the 
benefit to be had from a NNLL resummation. In order for the comparison to be consistent 
with HqT, which has just a single matching scheme, our uncertainty band here includes 
only scale variation with scheme (a), but not the central values from the other schemes. At 
low values of pt, veto ^ 20 GeV, it appears that there are clear benefits to be had from the 
NNLL resummation, with a significant reduction in the uncertainty as compared to a NLL 
result. On the other hand, for the practically relevant region, pt.vcto ~ 30 GeV, the bands 
are quite similar in size. This is perhaps not surprising since at this scale resummation only 
just starts to become relevant. 

The modest size of the differences between a transverse momentum veto on the Higgs- 
boson and one on the jet is to be expected given that the two observables differ only from 
O (a^) onwards. An interesting cross-check of the consistency between the two approaches 
can be obtained by taking the HqT calculation of the Higgs transverse-momentum and cor- 
recting it with the a^ difference between Higgs and jet transverse momentum distributions. 
Some care is needed in order not to introduce divergences at low pt, and we use the following 
matching procedure 

^(HqT+H2j@LO)/ \ - f^ , ^2,jet(Pt,veto) - S2,Higgs(Pt,veto) ^ ,(HqT)/ n /^ -, n, 

^jet [Pt,vcto) — [^~^ ) %iggs IPt.vctoj , I'-ij 

involving the difference of S2 terms for the jet-pt and Higgs-pt calculations normalised to 
aQ, as well as the Higgs pt-veto efficiency, CHiggs (Pt.vcto)- While such a matching does break 
the resummation accuracy, one may still expect it to have some meaning in the region 
of intermediate transverse momenta. The result is shown in Fig. 6 (right) as the band 
labelled "HqT matched with H2j@L0", using the MCFM calculation of H+2-jets at LO. 
It is in reasonable agreement with our NLL+NNLO calculation. This helps illustrate the 
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jet veto V. Higgs p, veto (HqT) 



jet veto: NLL+NNLO v. HqT with a^ correction 
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Figure 6. Left: comparison of our result for the jet veto efficiency with the calculation of a Higgs- 
boson transverse momentum veto efficiency obtained with HqT. For the purpose of the comparison, 
since HqT provides only a single matching scheme, we restrict the uncertainty band on our results 
to use just matching scheme (a). Right: comparison of our jet veto result with the jet veto result 
obtained by correcting HqT with the (relative) order a^ Higgs+2-jet contribution, as in Eq. (7.1). 
Again our band includes just scheme (a). 

consistency between the different tools. The study can be extended to NLO in the H+2- 
jet calculation [19], and the results appear broadly similar, though care is needed in the 
treatment of technical infrared cutoffs in the H+2-jet calculation. 

Let us now turn to POWHEG.^ Here we use default renormalisation and factorisation 
scales equal to M/f/2, as in our fixed order and resummed calculations. Fig. 7 (left) 
shows the results, showered with 3 different commonly used tunes of PYTHIA, Perugia 2011, 
Z2 [66] and AMBTl [67]. For the pt values of interest, all three are contained within the 
NLL+NNLO uncertainty band, though systematically above its central value. 

To probe the uncertainties in POWHEG plus PYTHIA, we followed the suggestion in 
refs. [10, 68], and varied renormalisation and factorisation scale independently around Mnf^ 
as described in Eq.(3.14) and fixed the parameter hfact to h = M/f/1.2. With these 
choices the Higgs transverse momentum spectrum of POWHEG+PYTHIA yields a reasonable 
agreement in shape with that of HqT. We recall that in POWHEG it is possible to split 
the real radiation contribution R into a shower part R^ and a finite one i?-^, with the 
requirement that R^ —^ R at small transverse momenta. Choosing a finite hfact (rather 
than hfact = oo) means that R'^ = h?/{h? + p^ ^)R and R^ = p\fj/{h? + plj. ^)R. 
Therefore this choice mainly affects the large transverse momentum region. Since the 
transverse momentum spectrum of shower events is not affected by scale variations, while 



^We used revision 1683. 
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NLL+NNLO V. POWHEG + Pythia 6.4 (3 tunes) 



NLL+NNLO V. HqT-rescaled POWHEG 
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Figure 7. Left: comparison of powheg + pythia 6.426 with our NLL+NNLO predictions; results 
are shown for three Pythia tunes. One of them, Perugia 2011, is displayed as a band, formed by 
the envelope of renormalisation and factorisation scale variations (Eq. (3.14)). Right: comparison 
of the HqT-rescaled POWHEG + pythia combination (Perugia 2011 tune) with our NLL+NNLO 
results. POWHEG and HqT are always used with MSTW2008NNLO PDFs. 



that of finite events is, this splitting provides one way of assessing of the uncertainties in 
POWHEG. 

In Fig. 7 (left) this is illustrated as a band for the case where POWHEG is interfaced to 
PYTHIA with the Perugia 2011 tune. This band is significantly narrower than our uncer- 
tainty band and one also observes that for low values of pt, veto it does not quite encompass 
the different tunes. We tend to believe that the POWHEG uncertainty band represents an un- 
derestimate of the true uncertainties. This suggests that it might be valuable to investigate 
other sources of uncertainty, for example variations of hf act. 

The right-hand plot of Fig. 7 shows the impact of reweighting events so as to ensure that 
the final Higgs pt distribution agrees with the output of the NNLL+NNLO HqT program. 
In this procedure, every (originally unweighted) event i from POWHEG+PYTHIA is assigned 
a weight Wi 



daH-^T 



Wi 



dpt,H 




^^POWHEG+Pythia 



dpt,H 



Pt.Hy 



(7.2) 



(i) 



where p^ ^ is the transverse momentum of the Higgs boson in event i. The result of the 
weighting for the main tune we consider, Perugia 2011, is to increase the efficiency a little, 
bringing it close the upper edge of our (full) band. The reweighted POWHEG+PYTHIA band 
is broader than that before reweighting, but still smaller than ours. This is natural given 
that HqT considers just a single matching scheme as compared to our multiple schemes. 
Note also that the reweighting of POWHEG + PYTHIA with HqT breaks the NNLO accuracy 
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of the HqT calculation, unlike the procedure of Eq. (7.1). 

Finally, we remark that while here we have considered just the POWHEG+PYTHIA 
results for NLO-matched parton showers, comparisons with MC@NLO [30] showered with 
HERWIG [69] are to be found in Ref. [10]. 

8 Conclusions 

In this article we have examined resummations and matching to fixed order calculations for 
jet vetoes in Higgs-Boson and Z-boson production. The NLL resummation has a particu- 
larly simple form, reducing to just a Sudakov form factor. Such a simplification had been 
observed in the past for the two-jet rate with the Cambridge algorithm in e"'"e~ collisions. 
Here the result holds for a broad class of jet algorithms, independently of the jet radius. 

To make phenomenological use of our result it was necessary to match the resummation 
with a fixed-order calculation, especially because experimentally relevant jet- veto scales are 
at the edge of the region where resummation starts to become relevant. One may express 
the result either as a jet-vetoed cross section or a jet- veto efficiency. We chose to concentrate 
on the latter for two reasons: firstly, the efficiency is closely related to the concept behind 
a Sudakov resummation; secondly, the perturbative structure of the efficiency provides 
useful handles when estimating uncertainties. The uncertainties on the efficiency should, 
we believe, be essentially uncorrelated from those on the total cross section, making it 
straightforward to also derive uncertainties on vetoed cross sections. 

For Z-boson production, the jet-veto efficiency turns out to be rather stable. For 
Higgs production, however, even if the NLL+NNLO matching provides some degree of 
stabilisation of the results relative to just NNLO, the final uncertainties are not negligible, 
cf. table 3. The problem can be traced back to ambiguities in how one formulates the 
fixed-order result for the efficiency and is closely connected with the poor convergence of 
the perturbative series for the total cross section. Essentially, considering that the known 
NNLO corrections to the total Higgs cross section are about 20%, one can expect that the 
unknown NNLO corrections to the 1-jet rate can also easily be of order 20% (and distinct 
from the uncertainty on the total cross section, associated with unknown NNNLO terms), 
which then inevitably propagates into the jet-veto efficiency. This issue is generally not 
addressed in existing calculations, with the exception of Ref. [54], which approaches it from 
a somewhat different point of view. From the perspective of resummations, such effects 
are N'^LL, and so appear to be beyond current technology. Steps towards NNLL accuracy, 
specifically terms related to the choice of jet definition, are discussed in the appendix. 

Given the NLL+NNLO calculation, there were various phenomenological questions that 
could be addressed. On one hand, it was useful to establish how directly the result could be 
applied to data. Effects of potential relevance include finite detector acceptance and non- 
perturbative corrections. Based on Monte Carlo simulations, both appear to be modest 
compared to other remaining uncertainties. We also compared our results to the predic- 
tions from tools used by the LHC experiments, such as the NNLL+NNLO Higgs-boson 
Pt resummations from the HqT program and the NLO parton-shower generator POWHEG 
(with PYTHIA), with or without rescaling corrections from the HqT program. We observe a 
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pattern of slightly higher efficiencies with these tools than in our NLL+NNLO calculation, 
but all remain consistent to within uncertainties. 
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A NNLL i?-dependent terms 

The CAESAR paper [33] was primarily concerned with NLL resummation. Nevertheless, it 
did provide certain formulae that are of use in a NNLL context, as part of the verification 
of the particular classes of terms that could be neglected in a NLL calculation. 

It is our understanding that these formulae provide key pieces that are missing in order 
to relate existing NNLL resummations for the Higgs or vector boson pt distribution to a 
NNLL resummation of the jet veto efficiency and specifically the dependence of the NNLL 
corrections on the jet algorithm and radius. 

While it is beyond the scope of this article to provide a full NNLL analysis of the jet 
veto efficiency, we believe it is useful, nevertheless, to make the NNLL i?-dependent results 
available. 

There are two sources of correction, related to the fact that the emission of two soft- 
collinear gluons can be separated into two parts: an independent-emission term and a 
correlated emission term, the latter being non-zero only when the two gluons are close in 
rapidity: 

M%{h,k2) = Mlih) Ml{k2) + M%{h, k2) , (A.l) 

where Ai'i{k) is the factorised matrix element for the emission of a single gluon and 
A^g„(fci, A;2) is the correlated-emission part of two-gluon emission matrix element in the 
limit where both are soft and collinear with respect to the beam direction. There is also a 
corresponding correlated-emission component for the production of a quark-antiquark pair, 
A^qq. We will use Ai"^ to denote Ai'i„ + 2A^qq. For double soft (and collinear) emission 
off a quark line, the non-trivial correlated part of the matrix element was given as early as 
Refs. [70, 71] and was rederived for more general event structures in Refs. [72, 73]. Apart 
from colour factors (C^^, CfCa, CpT^nf respectively for independent and correlated gluon 
and quark-pair emission from a quark line; C\, C\ and CATRUf for emission from a gluon 
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line), the matrix elements are the same regardless of whether the gluon pair or qq pair are 
emitted from a quark line or a gluon line, as can be seen clearly in the formulae of Ref. [73]. 

A.l Correlated emission 

We start with the correlated-emission component, and in particular from Eq. (D.8) of the 
CAESAR paper, which provides the NNLL contribution associated with the presence of any 
number of independent emissions and one correlated pair: 



J^' 



:orrel 



exp ( - / [dk]Mg^^i.{k) j x 

.^[\J\j idh]Mlrciki) I ^ J [dka] [dkb]Mlg^,^{ka, h) x 
X [Q{v - V{ki, ...,kn, ka, kb)) - Q{v - V{ki, ...,kn,ka + h))] . 



(A.2) 



This formula can be understood as the correction that arises when an observable is sensitive 
to the kinematics of the individual a and b partons rather than just their sum. For the 
resummation of a sufficiently inclusive quantity, notably the boson pt distribution, it is 
explicitly zero. Consequently, the calculation of this term should allow one to relate the 
jet-veto resummation to the boson pt resummation.^ 

For compactness we have not explicitly written the {p} argument to V{. . .) in Eq. (A.2). 
We work in a limit where a^ <C 1, In 1/v S> 1 and a^ In 1/v is finite. The parameter e, which 
serves as a regularisation cutoff, is to be taken e <C 1, but also such that a^lnl/e ^1. In 
this limit the phase space integrals [dki] essentially extend up to a rapidity \yi\ < In 1/v and 
the dependence of the precise upper rapidity limit on a given parton's transverse momentum 
will turn out to be irrelevant at our accuracy. Additionally, to within a factor 0{e), all 
emissions have a pt ~ vMb = Pt,veto (v = Pt,veto/MB and Mb, we recall, is the boson mass). 
The matrix elements include a subscript ^^rc" to indicate that the strong coupling is to be 
evaluated at a scale of order vMb ■ 

A NNLL contribution arises only when all emissions i = 1 . . .n are well separated from 
each other and from a and b. For the case of the hardest jet's pt, we then have 



Q{v - V{ki, . . . , kn, ka, kb)) 

e{v-V{ki,...,kn,ka + kb)) 
Making use of the fact that 



n 

[[e{v- 


-v{h)) 


n 

[[e{v- 


-v{h)) 



Q{v-V{kaM), (A.3a) 

Q{v - V{ka + h)) . (A.3b) 



E ^ ( n / [dk^^Ml,,{k)Q{v - V{h))] = exp ( T 

n=0 Vi=l-^™ / ^-^^ 



[dk]Ml,M 



(A.4) 



^We will not discuss here the effects of hard coUinear emission off two different legs [74], which can 
introduce an extra source of differences between boson and jet pt resuininations. 
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we can then rewrite Eq. (A. 2) simply as 

_^correl ^ ^ J [dka][dh]M%^rcika, h) [Q{v - V{ka, h)) - Q{v - V{ka + kb))] , (A.5) 

i.e. it reduces to a pure two-gluon result (with a running-coupling). 

It is simplest to first evaluate the leading /^-dependence of this formula in the limit of 
a small jet radius R. When the angle between the two partons a and b is small, A^^ = 
iUa — Vb)^ + (</'a — 'Pb)'^ ^ Ij we Can write the phase-space and matrix element as 



2Cas{kt^ab) dkt^ab d4>ab 

t,ab 
as{h,ab^ab) dAab 



[dka\[dkb\Mgg^^^{ka,kb) = — dyab^ 



dz{Pgg{z) + 2nfPgg{z)) , (A.6) 



where C = Cp or Ca depending on the nature of the beam partons, kab = ka + kb and the 
Pgg and Pqg are the real parts of the usual leading order splitting functions: 

Pgg{z) = 2Ca (-^^ + i^ + z{l - z)\ , (A.7a) 

P,g{z) = Tr [z^ + (1 - zf) . (A.7b) 

With these variables, for A^t *C 1, the difference of functions in Eq. (A.5) reduces to 

[e(^ - v{ka, h)) - e{v - v{ka + h))] = 

Q{Aab - R) [Q{vMb - max(z, 1 - z)kt^ab) - Q{vMb - h^ab)] , (A.8) 

as long as one restricts one's attention to jet algorithms from the generalised- /cj family. We 
can now rewrite Eq. (A.5) as 

X ^Ah^h^dZ {Pgg{z) + 2nfP,g{z)) X 



X 



(^hft [Q(yMB - max(z, 1 - z)kt,ab) - Q{vMb - kt,ab)]) ■ (A.9) 



A first point is that kt^ab will be limited to be of order vMb = Pt,veto- Secondly, while 
we take R < A^t <^ 1, we will still assume that a^lni? is negligible. Accordingly we can 
replace each of the running couplings with as{pt,vcto)- This puts us in a position to carry 
out the integrations in each of the three lines of Eq. (A.9) independently. The contents of 
the first round brackets on the first line give R'^ = 4Cas (pt,vcto )/'?'" In 1/v; the second set 
of round brackets on that line gives — Ini? (for now we neglect the O (1) contribution from 
the ill-defined upper limit in Aab)', the last line gives ln(l/max(z, 1 — z)). We are therefore 
left with 

^— -R'S(14^0 (1)) -%^ I' .4 iPggiz) + 2n,P,giz)) lu --^ . 

(A.IO) 
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This is straightforward to evaluate and gives 

^orrei „/ «s(pt.vcto) [^ 127r2 + 132 111 2 - 131 23 - 24 111 2 \ / 1 

^ = i?i. ^ [Ca ^^ + n^Tn j (^hi j^ + 0{l^ 

(A. 11a) 
~ R'^^£^^hl^ (1.09626C7A + 0.l76791n fTa) (in -^ + O {!)] . (A.llb) 

For a complete evaluation of the J^'^°'^^^^^ we start again from Eq. (A. 5). We multiply 
and divide A4'^gj.^{ka,kb) by a factor A4'^^^^{ka)M'^rc(.^b), and make use of 

[*lMj,„(^ = ?i^*£|l|*rf,. (A.12) 

We then replace the integration measure dkt[,/kt^hd4'bdyb with d(^/(^dA(j)dAy where ^ = 
kt^b/h,a, A(j) = (pb- (pa, Ay = yb- ya- This gives us 

-correl _ 1 f f^^sC dcPa ^^, \ 2a,C [^ dC [^ dA(b 



M'iJka-.kh) dki n 



Mjika)Mjikb) kt,a 

As before, the first factor in round brackets on the first line will integrate to give R'^. The 
difference of functions on the second line will be non-zero only when the two partons are 
separated by Acp"^ + Ay^ > R^ (this statement holds only for a recombination scheme such 
as the £^-scheme, which directly sums 4- vectors). In the ratio of matrix elements we drop 
the "re" subscript, since at our accuracy running coupling effects are entirely accounted 
for in the scale choice that we will make, asiPt,veto), for the explicit factors of as (cf. the 
discussion in the small- i? limit). Observing that kl^^/k'^^^ is independent of fef ^ and equal to 
{1 + C,"^ + 2C, cos A(/)), we can then perform the kt^a integration together with the 0-function 
constraints in the square brackets to yield a factor 



2 



l + C^-F2CcosA0' 
max{C^, 1} 



(A.14) 



Our result for J^'^°"'^^ is then 



TT Jo C J-n 27r y_oo 



1 Mlg{ka,kb) 1 
"" 2\Ml{ka)Ml{kb)2 '' 



1 + C + 2C cos A(j) 
max{C^, 1} 



(A.15) 



In both R'q and the explicit factor of a^, the coupling is evaluated at scale Pt,voto) as before. 
The derivative with respect to R of Eq. (A.15) can be straightforwardly evaluated as an 
expansion in powers of R. Integrating that expansion gives us a result for J^'=°'''''^l that 
is missing a constant of integration. This constant can easily be determined through a 
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, coefficient of R' C/^aj/Tc 



f'^°"^[ coefficient of R' n,ag/ji 



numerical result 
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Figure 8. Comparison of the numerical and small- i? analytical determinations of the Ca and Cp 

pieces of J^corrcl 

numerical integration of Eq. (A. 15), which also allows for a check of the range of validity 
of the power series in R. The result that we obtain is 



J^' 



;orrel 



R 



/ as{p 



't, veto J 



B 



vr 



-131 + 127r2 + 1321n2)C^ 1.74 (23-241n2)nf 0.84 

In ^^ + -^ — -^ In ■ 



72 



R 



72 



R 



+ 



(1429 + SeOOvr^ + 12480 In 2) Ca + (3071 - 16801n2)n/ ^ 
~^ 172800 

-9383279 - 11760071^ + 1972320 In 2) Ca + 2(178080 In 2 - 168401)n/ 



406425600 



i?^ 



+ 



(74801417 - 33384960 In 2)Ca + (7001023 - 5322240 In 2)n/ 



97542144000 
This result is compared to the full numerical determination in Fig. 8. 



R'^ + O {R^) 



(A.16) 



A. 2 Independent emission 

In the case of the jet observables that we have been discussing, the starting point for the 
independent-emission correction is 

jrindep ^ g^p f r [dk]Ml,^{k)\ X 

ZMn / \dh]^l,,{h)Q{v -V{h))\ 1 [[dka][dh]Ml,,{K)Ml,,{h)x 

n * \^ — 1 J ^V I'd 



n=0 



\i=l' 



X [e{v - V{ka, kb)) - 9(7; - Vika))e{v - V{kb))] . (A.17) 



The second term in square brackets on the last line corresponds to the approximation made 
in obtaining J-" = 1 at NLL; the first term corresponds to the actual value of the observable. 
The evaluation of this formula largely follows the working given above in the correlated 
emission case, with the main difference that now the only region that contributes is that 
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where the two gluons have A^;, < R. The result for i? < vr is 

_^mdep ^ _j^, as{pt,.eto)C 1 f""^^^^^^^ /"" ^ T ^ l^ f 1 + C^ + 2C COs(A,fc COS V^) 

^ 71" 2! io "^ "Vo 27r io C |_ maxjC^l} 

(A. 18a) 
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